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Abstract. 



The black hole model with a self-gravitating charged spherical symmetric 
dust thin shell as a source is considered. The Schroedinger-type equation for 
such a model is derived. This equation appeared to be a finite differences 
equation. A theory of such an equation is developed and general solution 
is found and investigated in details. The discrete spectrum of the bound 
state energy levels is obtained. All the eigenvalues appeared to be infinitely 
degenerate. The ground state wave functions are evaluated explicitly. The 
quantum black hole states are selected and investigated. It is shown that 
the obtained black hole mass spectrum is compatible with the existence of 
Hawking's radiation in the limit of low temperatures both for large and nearly 
extreme Reissner-Nordstrom black holes. The above mentioned infinite de- 
generacy of the mass (energy) eigenvalues may appeared helpful in resolving 
the well known information paradox in the black hole physics. 
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I. Introduction. 



The fate of black holes became a subject of interest since Hawking' s dis- 
covery that the large black holes should emit a blackbody radiation with 
the temperature depending on their parameters. The problem is that the 
temperature of the (uncharged, nonrotating) black hole grows when its mass 
decreases, and the density of the emitting radiation can no longer be con- 
sidered as negligible, thus requiring an account of the back reaction of the 
radiation on the black hole metric. To treat such a self-consistent problem 
is an extremely difficult task. 

Moreover, the values of densities and sizes of the objects involved are such 
that the classical gravity needs quantum gravity corrections. 

There is a general belief that the back reaction and/or the quantum grav- 
ity corrections would prevent the black holes from complete evaporation. 
This would lead also to numerous cosmological consequences. 

The full treatment requires the field theoretical or (super) string consid- 
eration. 

In the present paper we try a quite different approach. We consider a 
spherically symmetric black hole (both charged and uncharged) as a quantum 
mechanical object (like a hydrogen atom). We made the simplest possible 
choice for its source - a self-gravitating spherically symmetric dust shell. The 
thin shell limit in General Relativity elaborated in a geometrically invariant 
form by W. Israel [§] make many problems exactly solvable and, though very 
simple, allows one to obtain important physical results both in black hole 
physics |]14| and cosmology |jl5 . 



The thin shell model allows us to reduce the problem in the spherically 
symmetric case to the essentially the one dimensional problem which is rather 
easy to quantize. The quantization procedure used here was first proposed 
in the paper by V.A.Berezin, N.G.Kozimirov, V.A.Kuzmin and I.I.Tkachev 
[^] and then in |^. The discrete mass spectrum for the self-gravitating dust 
shells was suggested in 0. In a paper [|I^ it was generalized to electrically 
charged shells and used then to obtain a black hole spectrum, which appeared 
compatible with a Hawking's blackbody radiation in quasi-classical quasiclas- 
sical regime, i.e. for large as well as nearly extreme Reissner-Nordstrom black 
holes. 

In the present paper we summarize the results, obtained earlier and con- 
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struct for them a solid base. Namely, we describe a method of solving our 
Schroedinger equation in finite differences with Coulomb potential, found 
the general solution, solved the boundary conditions and proved the exis- 
tence of discrete spectrum for bound states, evaluate explicitly the ground 
state wave function and showed that all the mass (energy) eigenvalues are 
infinitely degenerate in the frozen formalism. 

The plan of the paper is the following. In Sect. II we present the classi- 
cal equations of motion of a self-gravitating charged dust shells in General 
Relativity and discuss some of their properties. This model is quantize in 
Sect. III. The resulting Schroedinger equation is not a differential equation 
but the equation in finite differences. We found the asymptotics of its so- 
lutions in singular points (near the origin and at the infinity) in Sect. IV. In 
Sect.V we present the general solution to the equation in the momentum 
representation and discuss its properties. A nontrivial transition to the co- 
ordinate representation is given in Sect. VI and the fundamental solution is 
introduced. This solution is investigated in details in Sect. VII. In Sect. VII 
the appropriate boundary conditions at the origin and infinity are written to 
ensure the Hamiltonian be self adjoint operator on the positive semi-axis. In 
Sect. IX the discrete spectrum for bound states is found, the new polynomi- 
als are introduced and investigated (in particular, the generating function for 
them is evaluated). The ground state wave function is presented in Sect.X. 
In Sect. XI the notion of quantum black holes is introduced and the black 
hole mass spectrum is derived. This spectrum is shown to be compatible 
with the existence of Hawking's radiation in Sect.XII. The Paper concludes 
with Discussions. 

Throughout the paper the system of units in which c = h = k = lis used, 
where c is a velocity of light, h is a. Planck's constant. A; is a Boltzmann's 
constant. 
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II. Classical Equations of Motion for a Self-gravitating 
Charged Spherically Symmetric Dust Shell. 

We start with a description of the model. Thus is just a self-gravitating 
spherically symmetric dust thin shell, endowed with a bare mass M and 
electric charge e. The whole spacetime is divided into three different regions: 
the inner part, the outer part and the 3-dimensional timelike hypersurface 
separating them. 

The general metric of a spherically symmetric spacetime has a form 

ds^ = Adt^ + 2Hdtdq + Bdq'^r^{t, q)dn^, (1) 

A > 0,5 < 0, 

where t and q are correspondingly timelike and spacelike coordinates. A, 
H and B are functions of t and q only, r{t,q) is the radius of a two-dimensional 
sphere (in the sense that the area of the sphere is inr'^O, 

dn^ = de^ + sin^ edifi^ (2) 

being the line interval of the unit sphere. 

Contrast to the flat spacetime, the normal vector to the surface r = const 
may be not only spacelike but also timelike. In the first case 

A = g'^^r^^r^p < 0, (3) 

and the corresponding region is called i?-region (here Qa/B is a metric tensor, 
its inverse, r^a denotes the partial derivative with respect to the cor- 
responding coordinates, Greek indices run from to 3). In the flat case the 
i?-region occupies the whole spacetime. In the timelike case 

A > 0. (4) 

Such a region is called the T-region (the notions of R- and T-regions were 
introduced in |[^). It easy to show that the condition r = ("dot" denotes 
the time derivative) cannot be satisfied in a T-region, hence it should be 
either r > (this region of inevitable expansion is called T+-region), or r < 
(inevitable contraction, a T_-region). Correspondingly, it is impossible to 
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get r' = ("prime" denotes the spatial derivative) in i?-regions, and a region 
with r' > is called an i?_-region, while that with r' < is an region. 
The i?+ and i?_ regions correspond to different sides of an Einstein-Rosen 
bridge (see, e.g., and references cited therein). 

The solution to Einstein equations representing the Reissner-Nordstrom 
(spherically symmetric, charged) black hole is well known and can be put in 
the form 

ds^ = fde - f-^dr^ - r^dn\ (5) 

where 

/ = 1 + — , 6 

and m is the total mass (energy) of the system, e is its electrical charge, k 
is the gravitational constant (equal to the inverse square of Planckian mass, 
K = M~i^, in the chosen units; note also that in these units the electric 
charge is dimensionless, and radius has dimension of inverse mass). For the 
zero charge case the above metric is reduced to the Schwarzschild solution. 
The Reissner-Nordstrom metric has been extensively studied (see for the 
detailed description and references). What is important for us is the follow- 
ing. 

The invariant A, introduced above, Eqn.(|^), is 

A = -f^-l + ^J^-^. (7) 

ry 

This function has two different positive roots, if ^/Km > |e|. 



= Km + V K^m^ — kc'. 



Km 



— vKprnF—Ke^ . (8) 



The surfaces r± = const are null surfaces separating regions with A < 
(i?-regions, r > r+, < r_ ) from that with A > (T- regions, r_ < r < r_(.) 
The r+ = const surface is the event horizon, the so-called outer horizon, like 
the event horizon in the Schwarzschild solution to which it reduces when e = 
0. The inner horizon, r_ = const, was recognized as a Cauchy horizon beyond 
which the solution cannot be continued unambiguously. There are infinitely 
many R± and T± regions. They form some kind of ladder on the Carter- 
Penrose diagram for a maximally analytically extended Reissner-Nordstrom 
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spacetime. But due to the generic instability of the Cauchy horizons only 
one portion (in between two nearest Cauchy horizons) can be considered as 
physical (for the details and references see and ). In many aspects this 
part is the same as the maximal analytical extension of the Schwarzschild 
metric. The main point for us is that there are two asymptotically flat R- 
regions (where r > r+), one of them, is the region external to the black 
hole, it is "our" side of the Einstein- Rosen bridge, the "other" side being 
i?_-region and is called also a wormhole region [0. 

The Reissner-Nordstrom black hole with y^m = |e| is called the extreme 
black hole. For ^/nm < \e\ no black hole solution exists. 

The metric of the 3-dimensional surface, S, representing the evolution of 
the thin shell can be written as 

ds^j. = dr^ - p\T)dQ\ (9) 

where r is the proper time measured by an observer comoving to the shell. 
This surface may be considered as embedded into 4-dimensional spacetime, 
then the Einstein equations on the shell may be expressed in terms of the 
extrinsic curvature tensor Kij, the curvature of embedding (Latin indices 
take values 0,2 and 3), and surface energy tensor Sij 0: 

[Ki]-5i[Kl] = -8nKSi, (10) 

where 

[Kn = Kiio^d)-Ki{^n) (11) 

is the jump in the extrinsic curvature across the shell (the coordinate, n, 
normal to the surface is chosen increasing in the outward direction). The 
remaining Einstein equations, (^^^ and " ^ components, may be written in 
the form (vertical bar denotes covariant differentiation with respect to the 
metric on S): 

siy + [Tn = 0, (12) 

K^Si + = 0, (13) 

K} = ^{K}{out) + K^M) (14) 

and coincide with the covariant energy-momentum conservation equation 
written for the shell ||^. In the case of the spherical symmetry Eqns.(|IO|-|l^) 
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can be expressed in terms of invariants and take the form ^ (note that now 
A2 = A3 and Aq = inv for any 3-tensor on E ): 

[Ki] = AttkSI (15) 
K] + [K^2] = 8^>^Sl (16) 



where 



4° + 2^(5° -5,2) + [To-] = 0, (17) 
i^o°--^(p+^-4..pTA (18) 



K', = -^^p^-A. (19) 

Here "dot" denotes the derivative with respect to the proper time r, cr = ±1 
(its meaning will be explained in a moment), and and are normal 
flow and normal stress components of the energy-momentum tensor external 
(internal) to the shell. For the Coulomb field 

T° = Tl = T,' = (20) 



for any choice of the (spherical) coordinate system. For a dust shell 5*1 = 0, 
by definition. So, 

- ^ (21) 

where M is identified as the bare mass of the shell. As can be easily shown, 
Eqn ([T6|) is a differential consequence of the Eqns.(|^) and (|l^, and we get 



eventually the single equation (inserting for A Eqn.(|^ and remembering that 
inside the shell the spacetime is empty) 



/— 7 /.„ ^ 2Km Ke^ kM 

a,Jp^ + l-aout\ P^ + l + ^ = 22 

V V p p^ p 

which is the main subject of our future consideration. 

The quantity a appeared for the first time in the expression for extrinsic 
curvature tensor, Eqns.(0) and (|19|). Its meaning is the following, a = +1 



7 



if radii increase in the outward normal direction to the shell, and cr = —1 if 
radii decrease. Therefore, in the i?-regions a does not change its sign, the 
latter being the sign of an i?-region. It is clear now that on "our" side of 
the Einstein- Rosen bridge we have cr = +1, this we shall call the "black hole 
case", while on the "other" side a = — 1, this we shall call the "wormhole 
case". In the flat (Minkowski" spacetime we have everywhere a = +1. 
Solving Eqn.(|2^) for m we get 

2 ^2 



m = M^^l- ^^^^ ' , (23) 

and now it is seen that this is nothing more but the energy conservation 
equation, the square root being the Lorentz factor written in terms of proper 
time derivatives. 

We shall be interested in the bound motion only. Let po denotes the value 
of radius of the shell at the moment of rest and po > r+ (i-e. we are outside 
the event horizon). Then it follows from Eqns. (|2^) and (p3D that 



m = M , (24) 

2po 

/ kM\ , . 

(^out = sign 1 , (25) 

1^ = 1-^. (26) 

It is seen now that the total mass m increases with bare mass M when 
<^out = +1, i-e. in the black hole case, and it increases when aout = —1 (the 
wormhole case). In the black hole case 

ml , . 

thus ranging from 1/2 (M — oo ) to 1 (M e/y^, and m ej as well 
). The opposite inequality is valid in the wormhole case. 

To complete our discussion let us consider the thin shell with negative 
total mass, m < o (but positive bare mass M), From Eqn. (^2]) it is readily 
seen that aout = — 1 everywhere, and the spacetime topology is 5*^ x rather 
than x R^ as it is for the positive mass m. Due to different topologies the 
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cases m > and m < in quantum theory should be treated separately. In 
what follows we will consider only positive masses, m < 0. 



9 



III. Quantization. 



To quantize a system we have to construct first the classical Hamiltonian 
and then to replace the classical Poisson brackets by the quantum commu- 
tation relations. In making all these we will follow the procedure elaborated 
in § and §. 

Since m is a total mass of the system and Eqn.(^) represents the en- 
ergy conservation law it is reasonable to consider it numerically equal to the 
Hamiltonian of the system in question. Then the conjugate momentum p, 
Lagrangian L and Hamiltonian H can be found from the following relations 

dL , , 

31 

H = pp-L = p—-L, (29) 
which give for the conjugate momentum 
Substituting for H the Eqn.(p3|) one gets 



p = M\og{p+Jp^ + l) + F{p), (32) 



L = M(plog(p + - + PF^P)^ (33) 

where F{p) is an arbitrary function. The choice of this function does not af- 
fect the Lagrangian equations of motion, and leads to the canonically equiv- 
alent systems for Hamiltonian approach. This choice may become important 
when one explores the rather complex structure of maximal analytical exten- 
sions of black hole spacetimes trying to construct a wave function covering 
the whole manifold as was shown for a specific example in 0. But here we 
do not need such a complication and F{p) = 0. From Eqn.(p2D we get for p 
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P = sinh^, (34) 
and after inserting this into Eqn.(p3D, the Hamiltonian H is 

i7 = McoshA_^i^^!^. (35) 

It is convenient to make a canonical transformation to the new variables (M 
is the bare mass of the shell) 

x = Mp, n = ^p, (36) 

then 

H = M (^cosh n - • (37) 

Let us now impose the quantum commutation relation 

[n, x] = -i. (38) 
Then, using the coordinate representation with 11 = —id/dx and the relation 

e-'^m{x) = m{x -i) (39) 
we obtain the Schroedinger stationary equation 

H^{x) = E^{x), (40) 

— ( {"^{x + i) + "^{x - z) \E'(x) j = m'^{x). (41) 

Introducing the notation e = m/M we can rewrite this as 

^(x + z) + ^(a;-i) = ( 2e+ j ^(x). (42) 

Some remarks are in order. First, the Schroedinger equation obtained is 
not the differential, but the equation in finite differences. This is a direct 
consequence of the quantization in proper time which is quite natural in the 
framework of General Relativity. Second, the total mass m is not arbitrary 
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anymore but it is an eigenvalue of Hamiltonian operator, subject to the 
condition that for corresponding eigenfunctions the Hamiltonian is Hermitian 
on the positive semi-axes. Third, the Hamiltonian picture we got may not be 
equivalent to that which could be obtained directly from the Einstein-Hilbert 
action because of making use of the proper time definition of which depends 
on the history of the shell. It would definitely inequivalent if we would 
start from the Einstein-Hilbert action and express it in terms of proper time 
before making variations. But we made it the other way around. We started 
from equations of motion written already in proper time and constructed 
the new Lagrangian which results in just the equations in proper time with 
no restrictions or/and constraints. The procedure is thus self-consistent. 
We believe then that the resulting mass spectrum will be correct because it 
depends on the existence or nonexistence of the appropriate solutions but 
not on their detailed behavior. 

It should be noted also that the step in our finite difference equation is 
along the imaginary axes, so the solutions should have certain analyticity to 
be continued to the complex plane. 

With all dimensional quantities restored the shifted argument {x ± i) 
becomes 

x±^^M(p±^^) ^p±—z, (43) 



M J Mc 

so for p ^ h/Mc we obtain the nonrelativistic limit. Expanding Eqn. (^2|) up 
to the second order in {Ti/Mc) we get 

1 <f^ KM^-e\ 

-m = {E- M)m (44) 



2M rfp2 2p 

which is just the nonrelativistic Schroedinger equation for radial wave func- 
tion in s-wave, {E — M) being the nonrelativistic energy of the system. For 
negative values of {E — M) acceptable solutions exist only for a discrete 
spectrum of energies, namely 

{E-M)n = ^, n=l,2,... 45 

For K = it reduces to the well known Rydberg formula for the hydrogen 
atom. 
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IV. Asymptotics. 



In this Section we will find the asymptotics for the solution of the finite 
differences Schroedinger equation in coordinate representation, i.e. for the 
Eqn.(ll). 

Our strategy will be the following. Expand, first, the left had side of 
Eqn.(^) into series in derivatives. We get 

■2 '3 

^(x + i) = ^(x) + t^' + + + ■■■ (46) 

"2 3 

^(x -i) = ^(x) - i^' + - ^^(3) ^ . . . ^47^ 
We truncate this series to obtain 



K-- 



= - E (-l)"^-^!^*^^^^ + (-l)^(2iV)! - 1 + ^) ^(^)' (50) 

with N arbitrary and a = kM'^ — . 

The Eqn.(|SD|) will be our starting point. To find the asymptotic expansion 
of the solution to this equation we shall use the matrix method described in 
the book First of all, we should transform Eqn.(^) into a matrix form. 
This can be made by redefinition 

^(^^^ = y'2N = -E^=l(-l)^-^gfTl/2^+i + (-l)^(2iV)! (e - 1 + f ) y,. 

(51) 

Forming from ?/j the one-row matrix we get 

Y' = Aix)Y, Aix) = Ao + -Ai, (52) 

00 
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where 



Y 



An 



0'2K+1 



V y2N J 
I 



1 
1 





\ fll 

l)^(2iV)!(e- 1) 



\N-K- 



.i(2A0! 
( 



a 





V 1 



\ 





1 

'^2N ) 



a2K = 0; 



K 



l,...,iV-l; 



\ 









/ 



(53) 



The Eqn.(^) has two singular points. The first one, at x = 0, is a regular 
singular point (for all the definitions and theorem proofs see ) . The second 
one, at X = oo, is an irregular singular point. Their treatments need separate 
considerations. 

Let us start with a regular singular point at x = 0. The matrix equation 
is usually presented in the following form in this case 



xY' = B(x)Y, 



(54) 



B{x) =Bo + xBi, 



(55) 



-Bo = Ai] Bi 



An 



(56) 



Because the matrix Bq has only zero eigenvalues it can not be transformed to 
the diagonal form, but only to the Jordan's normal form by some similarity 
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transformation. The matrix T corresponding to such a transformation is not 
unique but this does not affect the result. The Jordan's normal form of the 
transformed matrix is unique. Thus, with Y — TZ, we have 

xZ'^CZ; C^Cq + xCi] 
Co = T-'AoT = H; 



T 



( 




E 



/ 1 




1 
1 



V 7 



\ 




V 



\ 




/ 
/ 



/ 





1 


. 


. 










1 . 


. 







I 


. 


. 1 


/ 






C2 


C3 


C4 . 










1 


. 


.. 











1 . 


.. 











. 


. . 1 


7 








. 


.. 



_ 2{-l)^+i 
C2ft:+2 - (2A')!a ' 



7 = (-l)^(27V)!|. 

Now the solution can be written as follows 

y = T ( / + £ x'-pJ (/ + //logx), 



r=l 



(57) 
(58) 



where J is the unit matrix, and the matrices can be calculated from the 
following recurrent equation 



UPr - PrH - rPr = -CiP, 



r—li 



(59) 
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Po = I, r = l,2, 



The matrix Y is to be understood in that sense that each of its row gives a 
hnearly independent fundamental solution to the Eqn.(|52D. It is not difficult 
to see the logarithm will appear in the first line (which is = \1/ and the 
only one of interest) starting from (1 + log a;)-term, and when N —>■ oo 
it disappears leaving 1 (this is an unexpected and remarkable result). 

The final result is that the wave function has the following asymptotics 
at X — (we are writing down only the leading terms) 



n = 0,1,2,.... 



(60) 



Let us now proceed with the asymptotical behavior of the solution near 
the irregular singular point X = oo. The procedure will be different in this 
case. We start again from Eqn.(^) and make a similarity transformation. 
But now the matrix Aq plays the role of -Bo(= ^i), and it is possible to trans- 
form it to the diagonal form (assuming that all its eigenvalues are different 
what is, fortunately, the case as will be clear in a moment). With Y = TZ 
we have now 



Z' = cz- 



C — Cq -\ — Ci, 

X 



Co = T MoT = diag(Ai, A2Ar); 



T 



/ 1 

Ai 



1 

A2 



2N-1 _)^2Af-l 



1 \ 

^2N 



.. A 



2N-1 



(61) 



Denoting by tij the elements of the inverse matrix T ^ we obtain for Ci. 



C, = {-lf{2N) 



t 



1,2N 



t2N,2N 

SO the matrix Ci has all the rows identical. 




(62) 
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The characteristic equation for A reads as follows 



N _^Ky2K 

which in the hmit N ^ oo becomes 

TTl 

cos A — e— — , (64) 

and it is clear now that all the eigenvalues are different. Note also that for 
e < 1 (we restrict ourselves to the positive values of mass) the eigenvalues 
are real, while for e > 1 they are pure imaginary. 
After one more transformation 

Z = P{x)W, 
W' = {P-^CP - P-^P')W = DW, 

oo 

r=l 

oo 

D = Co + Y. ^rX-\ (65) 



where all matrices can be chosen to be diagonal, and the diagonal elements 
of matrices P^ are all zero, the solution can written as follows 

Y = TP{x) exp |y (^Co + ^ DrX-'^ rfxj 

= r(7 + F(x))expCox + Dilogx, (66) 

where the matrix F[x) is holomorphic at X = oo. All the matrices Pr, 
and, hence, F can be calculated one by one from a recurrent equation readily 
derived from the above written expressions; they are of no importance for 
us. We are interested only in the matrix Di the equation for which reads as 
follows 

= Ci + CoPi - PiCo, (67) 
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and we are looking for the solution with Di having only diagonal elements 
nonzero and Pi having zero diagonal elements. It is easy to see that the Di 
consists of diagonal elements if the matrix Ci, hence, 

Di = {-lf{2Ny.^dmg{ii,2N,i2,2N, W^), (68) 

so the leading term in the first line of the matrix Y (remember we need only 
these elements) is 

^ gA.Xg(-l)'^(2JV)!(a/2)?,,2ivlogx^ j^gg^ 
S = 1,2, ...,2N 

where s is numerating linearly independent fundamental solutions. 

The only thing left is to calculate the second exponent in Eqn. (|70D and 
take a limit N ^ oo. Without loss of generality we can do this for only one 
value of s, let it be s = 1. Making use of a specific structure of similarity 
matrix T, Eqn.(pT|), we get for ii^2N 



'1.2N 



(70) 



(Ai-A2)(Ai-A3)...(Ai-A2^)' 
But from Eqn.(^) it follows that 

R2N = — {2K)\ ^ ^ {2N)\ ~ ^^^^^ - A2) ... (A - A2Ar), 

i^. = ^''^^^'^ = ' ~ ^^^^ 

Now it easy to take the limit ^ 00. We get finally 

^ = e^^e-*'°sx, (72) 

cos A = e. 

For e > 1 the eigenvalues are pure imaginary, A = ±i | A | , and we have ingoing 
and outgoing waves at infinity. For e < 1 the eigenvalues are real, 

A = ±cos"^ 6 + 2^77, A; = 0,±1,±2,... (73) 
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Note that our original equation in finite differences has the following inter- 
esting property. Given the solution, say \l/o, we can construct new solution, 
multiplying by any function C{x) which is periodical with the period 
equal to the imaginary unit, i. That is, 

^0(2^) - solution 
*i(a;) = C{x)'^o{x) - solution, if 

C{x + i) = C{x). (74) 
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V. Solution of the Schroedinger Equation in Momen- 
tum Representation. 



It is convenient to introduce the following parameters, A and 

e = cos A, a = kM^ - = 2(3 sin A, (75) 
in terms of which our Schroedinger equation becomes 

^{x + i) + ^{x-i) = 2 ^cosA + ^^^^ ^(a^)- (76) 

Since for the black holes a > 0, the signs of A and beta are the same. We 
choose A > and, hence, /? > 0. 

Let us return to the operator form of the equation, 

(^coshp — /5sin Ax~^^ = cos A\l/. (77) 

In the momentum representation operators p and x act as follows 

d , ^ 

p^'p = p^'p, x^p = (78) 

where|Psip is a wave function in the momentum representation. 

The operator x~^ is not well defined. The ambiguity can be removed by 
adding suitable terms to the potential which are proportional to 5-function 
and its derivatives at the origin. But, instead, we can multiply the equation 
by operator x from the left. By doing this we get 

d 

i— (coshp — cos A)\l/p = /3sin A\l/p, (79) 
op 

or 

d , ^ sinh p + i3 sin \ . , 

— log % = f ^— . 80 

op cosh p — cos A 

After introducing a new variable, z = e^, the Eqn.(^OD takes the form 

9 , ^ z^ + 2iPsm\zi 1 P + 1 (3-1 

T^log'^p = — — r — — T = + -, 

oz z{z^ — / COS A^; + 1) z z — zq z — zq 

where zq = e*^, Zq = e"*'^. (81) 
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The above equation can be easily solved, the result is 



Z ( Z — Zq'^^ 



\Z — Zq)\Z — Zq) \Z — Zq/ 



z = e^. 



This solution has very important property, it is periodical with the pure 
imaginary period 27ri. This property will be explored in the next Section. 
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VI. Transition to the Coordinate Representation. Fun- 
damental Solution. 



In this Section we will transform the solution found in the momentum rep- 
resentation, to the coordinate representation. It can be done by the inverse 
Fourier transform, 

oo 

^(x) = f e'P'^'^pdp. (83) 



^27r 

— oo 

In the previous Section we obtained the unique (up to the multiplicative 
constant) solution for \E'p , namely, Eqn.(|83D. But due to the periodicity of 
\I'p we can shift the argument p —>■ p + 27Tki{k = 0, ±1, ±2, ...) which will 
result in the shift of the path of integration in the complex momentum plane 
from the real axis to the parallel one. And after such a shift we again obtain 
a solution in the coordinate representation. But, 

oo 



^,{x) = ^ / e'P^%dp 

— CO 

oo 

vl/2(x) = -= / e'^P+'^'''^''^^p+2.k^dp = 

—oo 

\/2tt J 



27r 

— oo 

^2(0;) = e-'"'="^i(a;). (84) 

Thus, given one solution (say, \&o) ? we can construct in this way a count- 
able number of solutions solutions, and, in general, 

^generalix) = ( £ Cfce'^-^'^ | ^o{x). (85) 
\k=—oo I 

The infinite sum in parenthesis is nothing more but a Fourier series for a pe- 
riodical function with an imaginary unit period i. In this way we reproduced 
the property of the solutions in coordinate representation stated in Sect. 4. 
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How many different solutions do we need to construct the general so- 
lution? The very fact that in the momentum representation we obtained 
essentially one solution proves that for this we need only one solution 
which we call the fundamental solution (though it deserves the name super- 
fundamental). In the rest of this chapter we will construct the particular 
fundamental solution by the suitable choice of the contour of integration in 
the inverse Fourier transform integral. 

Note, first of all, that our solution, Eqn.(|83D, in momentum representation 
has countable number on branching points in the complex momentum plane 
which can be combined in pairs, {iX+2TTki, —iX+2TT{k+l)i), k = 0, ±1, ±2, .... 
Connecting two branching points in each pair by a cut we obtain the complex 
plane with countable number of cuts. On the corresponding Riemann surface 
our solution is a single valued analytical function of complex variable. 

Our choice of the contour of integration is as follows. We will integrate 
first along the real axes from left to right (i.e., from — oo to +oo), then 
along a short curve at the right infinity (p — > p + 27Ti,p — > +cxd), then along 
the straight line y = 27ri, parallel to the real axis, from left to right, and 
finally along a short curve at the left infinity back to the real axis. The 
integration along the short curves at infinities between the straight lines 
gives zero contributions for positive values of x (for negative x we can choose 
the straight line y = — 27ri instead of y = 2tti). The integration along each 
straight line gives us a solution the linear combination of which is again a 
solution. Thus, the inverse Fourier integral along such a closed contour gives 
us a solution. This contour can be distorted to become a contour around the 
cut (Ai, {2ti — for X > (or around the cut (— Ai, — (27r — for x < 0). 



Thus 




(86) 




C- 



(87) 



where 




(88) 
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In what follows we restrict ourselves to the case of positive x only. 
The above integral representation for can be simplified if we integrate 
Eqn.(^) by parts. Using the fact that 

zdp _ 1 ^ z-zq 

{z — zo){z — zq) 2isinA V z — zq 

we get 

c+ ^ ^° 

(the extra term vanishes because our contour has no boundary and for con- 
venience we omitted some constant factors). 

We will use use Eqn.(^) as an integral representation for the fundamental 
solution of our finite differences equation, Eqn.(main). 
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VII. Investigation of the Fundamental Solution. 



Let us investigate the fundamental solution, Eqn.(pO|), in more details. 
Our aim is to reduce the integration along the closed contour C+ around the 
cut to the finite interval between the corresponding branching points. But for 
/3 > 1 (remember that we have chosen (3 positive) we have the nonintegrable 
singularity at the lower integration limit. To avoid this difficulty we need 
some recurrent relation for lowering the parameter (3. To find such a relation 
we integrate by parts Eqn. (|90|) in the following way. 




iix + i)zo i e*f(^+*) 

J \Z — Zn, 

c+ 

-[x- i{f3 - l)^p-i{x) - e~'^^^p.i{x + i)] (91) 



Finally, we get 

IX 



^p{x) = ^^fs-iix) + J— J {"^beta-iix) - e-'^^^_i(x + z)} . (92) 

From the above relation it is easy to derive the structure of \E'/3 for general 
values of /3 > 1. Let us take f3 = n + P, with /? < 1. Then, ^f/j is the sum of 
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two terms, which of them is the product of some polynomial of n-th degree 
and of or "ifuidet^ix + i). Thus, we can proceed assuming [3 <1. 
First of all, consider separately the case (3 = 1. We have 



c+ 



= - (47re"^^sinAj xe"^"^. (93) 

From this it follows that all for positive integer (3 = n has the form 

^„ = P„(a;)e-^^ (94) 

where Pn{x) are some polynomials of n-th degree. These polynomials were 
first introduced in . 

In the case < /3 < 1 the integrand in the Eqn. (pOD is integrable on both 
ends of the cut in the complex momentum plane and we are able to convert 
the integral along the closed contour (7+ into the integral along the finite 
interval. We get the following, 

v|/^(x) = j e'P-fpdp = (l - e^-^) x^p{x), (95) 
c+ 

where 

iX ^ ' 

Changing the variables, 

e? -e^^ = -(2zsinA)i/, (97) 

we get 

$/3 = (l - e-^^^) e-^^^/ (l - (l - e-^^^) yj-' dy. (98) 







Comparing this with the well known integral representation for the Gauss's 
hypergeometric function, F(a, 6;c;2;), , 

1 ) 

F{a, h- c- z) = — / ^"-^(1 - t)-''-^(l - tzY'^dt, 

B{b,c-b) J 
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1 

B{x, y) = J - t)y-^dt, (99) 



we see that the integral in Eqn. (p8|) is the Gauss's hypergeometric function 
F{a, b; c; z) with the following values of parameters, 

a=l-ix, 6 = 1-/?; c = 2; z = l- e~^*^, (100) 

and 

5(6, c - 6) = 5(1 - A 1 + /3) = (101) 

sm vrp 

Finally, for the fundamental solution ^^/^(x) we get 
^^(x) = (-47r/3e-^^ sin A) xe-^^'F (l - zx, 1 - /5; 2; 1 - e'^^^) . (102) 

The above expression was derived for < /3, 1 only. We shaw now that 
the same is valid for any value /5 > 0. 

Let us check first the relation for (3=1. Indeed, from Eqn( p.02| ) it follows 
that 

^1 = - (a-kc-'^ sin a) xe-^^'F(l - ix, 0; 2; 1 - e'^*^). (103) 

But F(a, 0;c;z) = 1 and we see that Eqn. ( p^) coincides with Eqn(p3|) de- 
rived directly from our integral representation. 

Second, we show now that the recurrent relations, Eqn.(p^), for lowering 
index are direct consequences of the properties of Gauss's hypergeometric 
function. Indeed, let us write the following two of the well known Gauss's 
recurrent relations between hypergeometric functions with different values of 
parameters , 

c(c — bz — a)F{a, b; c; z) — c(c — a)F{a — 1, b; c; z) 
+abz{l - z)F{a + 1, 6 + 1; c + 1; z) = 

c(c — az — b)F{a, b; c; z) — c{c — b)F{a, b — l;c; z) 
+abz{l - z)F{a + 1, 6 + 1; c + 1; z) = 0. 
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Subtracting one from another, we get 

{c-b)F{a, 6-1; c; z) = {a-b){l-z)F{a, b; c; z) + {c-a)F{a-l, b; c; z). (104) 

Substituting our values of parameters we obtain 

/5F(1 -ix,l -/?;2; 1 -e~2*^ = (105) 
{p-ix- l)e-^'^F{l - zx, 1 - (/? - 1); 2; 1 - e'^*^) 
+ (1 + ix)F{l - i{x -i),l-{l3- 1); 2; 1 - e""'^). 

And, in terms of we have 

^>,{x) = |vl/,_i(x) - ^ (vl/,_i(x) - e^^^>,.^{x - .)) | . (106) 
Using the original finite differences equation we can rewrite this in the form 

IOC / \ 

■^0{x) = ^/3-i(x) + j-^ (^/3-i(a;) - e-'^^^„i(x + i)) (107) 

which is just the recurrent relation, Eqn. (|92|) , derived directly from the inte- 
gral representation. 

And, finally, it can be easily shown that the original equation in finite dif- 
ferences. Eqn. (|76D , is a direct consequence of the following Gauss's recurrent 
relation , 

(2a—c—az+bz)F{a, 6; c; z) + (c— a)F(a— 1, 6; c; z)+a{z—l)F{a+l, b; c; z) = 

(108) 

At the end of this section we consider the asymptotics of the general 
solution at X — > and a; — *■ oo . 

From the hypergeometric expansion 

and the form of the general solution 

oo 



2nx 

k = ~OD 

^0 oc xe~^^'F(l -ix,l- (3]2;1- e-2iX) 
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it follows immediately that the asymptotics at x — > are 

x\ r = l,2.... (110) 

This is just what we obtained earlier in Sect. IV. 

To obtain the asymptotics at x — oo we make use of the following trans- 
formation relation for the hypergeometric function ||^ 

^/ , N r(c)r(c — a — 6) ^, , 

F a, h- c- z) = -(F a, b; a + b - c + 1;1 - z) + 

r(c — a)r(c — 6) 

(1 - ^r''-' ^^'^rf''^rfi~ '^ F{c-a,c-b;c-a-b+l;l-z), 
where r(. . .) is the Euler's function. We get, then, 
= (-47r/5e-^^('^+^) sin a) a; 

Writing the hypergeometric function in the form 
and using the relation 

T{z)Til-z) = -^ (113) 
smvr^; 

we can rewrite in the following way 

47r/?e-^^(^+i) sin A 



r(i - /?)r(i + /3) e^^-*^^ - e-^^+^'^Z^ 

I fro r(l + fc-ia;-/?)fc! 



a/.^A. r(l+^a; + /c)r(l + fc + /j) ^_,,,, 
fr r(l + fc + ia;-|-/3)A;! 
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Using then the asymptotic form of the Euler function, 

T{z) = z'-^V27T(f){-), z^oo (114) 
z 

we find the asymptotic behavior of our fundamental solution ^^(x), 

- .-A. A. ^ r(i + fc + /3) ) 



to ^! 



But 



oo 



= r(l - /3)e-M/?-i)2/?-ief (^-1) sin A'^-\ 
and we obtain finally 

= -27rz/3e-^^e-^ {^TT^^"""'^ " r(T- !?)' ^"'"'^} '^^^^ ' 
X — > oo (115) 

Prom this it follows that the asymptotic behavior of the general solution is 
exactly the same as was found earlier. 
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VIII. Boundary Conditions and Conserved Current. 



In the preceding sections we found the general solution to the equation 
in finite differences with Coulomb potential, 

^(x + + *(a; - z) = 2 |^cos A + ^ sin ^(x). (116) 

Now we should remember that this equation is actually the radial Schroedin- 
ger equation for zero angular momentum, and the coordinate x runs from 
zero to infinity, i.e., x takes values in the positive semi-axis only. Thus the 
Hamiltonian, 

should be the selfadjoint operator on the positive semi-axis rather than on 
the whole real dictated by the quantum mechanics postulates. In 

addition, the wave function should be square integrable on the semi-axis. The 
corresponding extension of the above Hamiltonian was found by P.Hajicek 
in 0. 

It appeared the wave function should obey the following boundary con- 
ditions at the origin, x = , 

^(2«)(o) = 0, ri = 0,l,.... (118) 

That is, the function itself and all its even derivatives should be zero at the 
origin. 

The appearance of infinite number of conditions is due to the infinite 
order of the operator. 

P.Hajicek found also a conserved probability current J{x) for the Schroe- 
dinger equation in finite differences, Eqn. (|116|) . It is given as usual by the 
equation 

J'{x) = i - ^if^*) (119) 

so that 

(^*^)- + J' = (120) 
Writing the Hamiltonian (|117|) as follows 

^--E^^^^^Hx)-^ (121) 
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we obtain for J(x) , 



oo / 1 \k k 

k=o 1=1 



P.Hajicek showed that the boundary condition (|118|) imphes that 



J(0) = 0. (123) 

The boundary condition ( p.l8| ) is a direct generahzation of the boundary 
condition in the case of nonrelativistic Schroedinger equation which is of sec- 
ond order. In this particular case the boundary conditions ( |118| ) are reduced 
to the single one, \l/(0) = , which is together with the square integrability 
condition ensures existence of a discrete spectrum for bound states. And, in 
general, let us remind the results of Sect. IV, where we considered the trun- 
cated Hamiltonian of 2A^th order. It was shown that the asymptotics at the 
origin are 

^ ~ X, x\ x^^'\ 1 + x^^-^ log a; (124) 
and at the infinity they are of the form 

^ ^ giAfcXgafclogx^ ^ _ _ ^^25) 

where are real for bound states. The general solution is defined up to a 
normalization factor and depends, actually, on {2N — 1) arbitrary constants 
which are to be determined using boundary conditions. The square integra- 
bility condition reduces the number of unknown constants to (A^ — 1). But 
for the truncated operator of 2A^th order we have N conditions at the origin. 
They can only be satisfied for a discrete spectrum of eigenvalues. 



Note that the last of the asymptotics ( |124| ) is excluded by the boundary 
condition at the origin. 

We can reverse the above consideration and start from the asymptotic 
behavior at the origin. To satisfy the boundary conditions there we have to 
exclude the last asymptotics ( p. 24] ) and all the even terms in the expansion up 
to the (2iV — 2)th order. By doing this can determine N of (2A^ — 1) unknown 
constants. Thus we are left with only (A^ — 1) constants for N boundary 
conditions at infinity. But the situation is not so simple in the case of an 
infinite order operator. First of all, how to get rid of the logarithmic term in 
asymptotics? The problem is that if we first go to the infinite order limit, then 



32 



the logarithmic term unconditionally disappears (|x| < 1!) leaving 1 as an 
asymptotics. But the latter can be easily reproduced by some infinite linear 
superposition of the remained asymptotics. Subtracting one from another 
we would obtain the solution which would satisfy formally all the boundary 
conditions at the origin still having enough arbitrary constants to satisfy 
boundary conditions at infinity. 

But if we first differentiate 2N times the asymptotic solution containing 
logarithmic term we will get (1/x) term. Thus, requiring that the "good" 
solution should be infinitely differentiable we can reach our goal. Moreover, 
we have to require the analyticity of the wave function at least on the real 
axis. Thus is because we implicitly have used the analyticity of the solutions 
in transition from the differential equation of infinite order to the finite dif- 
ferences equation. The situation is not good at the infinity either. We saw 
that in the case of our particular equation we can obtain a new solution mul- 
tiplying the known one by exp {—2Tikx) . Choosing the large enough value 
of k we are able to convert the "bad" at infinity solution to the "good" one. 

We will see in the next section that it is the analyticity condition which 
solves this problem. 
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IX. Discrete Spectrum. Polynomials and Generating 
Function. 



We have shown in Sect. VII that the fundamental solution 
^^(x) = (-47r/5e-^^ sin A) xe'^^F (l - ix,l - (3;2;1 - e'^'^) . (126) 
can be written in the form 



X 



T{ix + /3)e^^^e-^^ ^ ^_ ^ ^_ 



r(i + zx)r(i + /5) 



F(l — ix,l — l3;l — ix — /3; e 



and the asymptotic behavior at x — oo is the following 

[ r(i + p r(i - ) X 

X oo (128) 

From this it is clear that the fundamental solution is an analytic function 
only if 

(3 = n (129) 

(remember that we chose jS > 0). This leads to the discrete spectrum for the 
total mass m. Indeed,we have 

f] = —— =n, a = kM\ 
2 sm A 

m 

e = — = cos A, 
M 



4n2 



(130) 



For j3 = n the hypergeometric series in Eqn( |126D terminates and we are left 
with the following fundamental solution, 

^„(x) = P„(x)e-^^ (131) 
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where P„'s are some polynomials of order n. They can be calculated directly 
from the definition, 



P^[x) = i^-Anne-'^ sin A j a;e-^^F(l - zx, 1 - n; 2; 1 - e"^*^) . (132) 

We can also make use of the recurrent relation found earlier which for our 
polynomials Pn takes the form 



Pn+l{x) = Pn{x) + - {Pn{x) - e'^'^Pnix + z)| . (133) 

In this way we are able to calculate Pn{x) step by step, starting from Pi : 

Pi (a;) = (-47re-^^ sin a) x. (134) 



But there exists yet another way to derive polynomials ( |131| ) . To find it 



we start with the following differential relation for hypergeometric function 



F{a,b;c;z) O 



(^) {^''^'~'FMc;z)} = ^±p-z'-'F{a,b + k;c;z), (135) 



which in terms of polynomials P„ reads as follows, 



For /c = 1 it becomes (remember that z = exp — 2zA) 

Q f iXn 1 jA(n-l) 

OA [n(smA)" J (smA)"+^ 
Using the relation 



P^_^{x) = -e'^'^Pr.+iix) + 2e'^ (^cos A - ^ sin A j P„(a;) (138) 

which follows from the following Gauss's recurrent relation 

{c-b)F{a,b-l;c; z) = {2b-c+ {a-b)z)F{a,b; c; z) +b{z -l)F{a,b+l; c; z), 

(139) 
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we get 



■Pn-i{x) 



dX Vn(sinA)"y (sinA)"+i 

- 2 cot A , . PJx) H ———PJx). 

(sin A)" ^ ^ n (sinA)'^ ^ ^ 

After rescaling the polynomials 

Pn{x) = (-1)'^ (e^^ sin a)" Pn{x) (140) 
we arrive at the following relation 

ne-''^Pr.^,ix) = ^ {e-''^Pn{x)) (141) 
which finally results in 

e--P..«W = l(^)"(e-P,W) (142) 

Pi{x) = (47rsin2 A)x. 



Writing 



and noting that 



we get 



e-^^=f^^^y' (143) 
V cot A - i / ^ ^ 



/• x2n -2At 1 ^ /COtA + z\ 

(sinA^)a;e ^^"^ = ^ : (144) 

^ ^ 2 9cotA IcotA-W ^ ^ 



Prom this it follows the form of the generating function for our (once more 
rescaled) polynomials Iln{x) = {l/n)Pn{x) 



oo 



*(^, ^) = E , = E ^n{x)z- 

n=0 n=0 



9\-,. + cot A + A 

27re^^^ ^ 146 

V 2; + cot A - z / ^ ^ 
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X. Wave Function of the Ground State. 



In this section we give an example of solving the boundary conditions and 
find the wave function of the ground state (i.e., for n = 1). In the preceding 
section the fundamental solution for n = 1 was found to be 

= xe"^^ (147) 

(we omitted the irrelevant here constant factor). The general solution obey- 
ing the boundary condition at infinity (exponential falloff at x ^ oo) has the 
form 

oo 

^i,en = xe-^"Ecfce-'"'"' (148) 

k=0 

and the coefficients Ck are to be determined by solving boundary conditions 
at the origin, 

^(20(0) = 0, / = 0,1,... (149) 

Differentiating Eqn. (|148|) 21 times we get the following infinite set of linear 
equations, 

oo 

2lY,Ck{X + 27rkf-^ = 0, 1 = 0,1,... (150) 

k=0 

The determinant of this system is identically zero because the first line con- 
sists of zeros. 

We can truncate this system at some specific value of /, calculate all the 
determinants and minors and then, after appropriate renormalization (say, 
putting Co = 1), take the limit I oo. In our simplest case of the ground 
state it is possible to make all the calculations up to the very end with the 
following result 



fc/ f(A) • k = 0,l,... (151) 



Inserting this into Eqn. (|148D we can write the ground state wave function 
(up to the normalization factor) in a very simple and nice form, 

oo 
k=0 
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oo 




- + 1 



) 



( 




xe 



(cosh TTX) T 



(152) 



It is easy to see that after shifting the parameter A ^ A + 27r/, where / 
is a positive integer, we again obtain a solution satisfying all the boundary 
conditions (technically such a solution can be obtained if put zero first / co- 
efficients Cm, m = 0, ...(/ — 1) in the infinite set of Unear equations considered 
above and normalize to unity q). Moreover, the introduced earlier conserved 
current is identically zero for any superposition (with complex coefficients) 
of these functions. This means that each eigenvalue of a total mass (energy) 
is infinitely degenerate. Clearly, the origin of such a degeneracy lies in freez- 
ing all but one dynamical degrees of freedom. Our experience in quantum 
mechanics suggests that after restoring the frozen degrees of freedom this 
degeneracy disappears. As a result we expect that the "true" spectrum for 
thin shell masses will depend on more than one quantum number. Due to 
infinite degeneracy of the "frozen" spectrum the "true" one will, in fact, be 
quasicontinuous. And, as a consequence, the Hawking's evaporation spec- 
trum will be quasicontinuous as well. The existence of the ground state with 
nonzero mass means that the evaporation process should stop when the mass 
of a black hole approaches the minimum possible value O.dMpi) . And 
we may conjecture that the final state of the black hole evaporation will be 
characterized by the minimal value of mass but different wave functions of a 
ground state depending on details of mass distribution and evaporation pro- 
cess. This unexpected feature of the black hole mass spectrum (its infinite 
degeneracy) reminds the problem of so called hidden parameters in quantum 
mechanics and may be helpful in resolving the well known information para- 
dox in a black hole physics. And at the end of this section we give the result 
of calculation for the norm of the ground state wave function. The scalar 
product of any two ground state wave function with parameters Ai and A2 is 
given by 
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' 2-"^ j/'ti^ cf2,^l. (153) 

2tt ) 



47r5/2 r (I + V ' 27r 



where Ci^^y) is a Riemann's function of two variables. Thus, for a wave 
function with given value of A we get 



^-sW2-"Fpi)Cf2,A|. (154) 
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XI. Quantum Black Holes. 



In the preceding sections we have got a complete solution to the Schroe- 
dinger equation in finite differences describing a quantum mechanical be- 
havior of the self-gravitating electrically charged spherically symmetric dust 
shell. For the bound states we found the discrete spectrum of mass which 
coincides with Dirac spectrum if we put zero the so called the so called ra- 
dial quantum number. Similarly to the case of hydrogen atom, to which our 
spectrum reduces in the nonrelativistic limit, we may conclude that the shell 
does not collapse without radiation. The starting classical situation is how- 
ever different. The classical hydrogen atom collapses radiating continuously 
while the classical spherically symmetric shell collapses without radiation. 
But without such a collapse the event horizon and thus the black hole can 
not be formed at all. Moreover, if we calculate the mean value of radius of 
the shell using the wave function we will find that at least for large values 
of principal quantum number it is far away of classical event horizon, so the 
shell spends most of its "lifetime" outside the classical black hole. So, to get 
a black hole solution we need some new criterium and more detailed investi- 
gation of the spectrum is needed. Let us consider the dependence of a total 
mass m on a bare mass M , the other variables, e and n fixed. We see that 
there are two branches, an increasing one for 

(kM^ - e^) (3kM2 - e^) < An^ (155) 

and a decreasing one in the opposite case. Note that now 

TTli /2 / 6^ 

and this value is greater than the corresponding classical value. 

The increasing branch corresponds to the "black hole case" while the 
decreasing branch - to the "wormhole case" . Using a "quasiclassical" argu- 
mentation we can say that for the states obeying inequality ( |155| ) the expec- 
tation value of the radius lies outside the event horizon on the "our" side 
of the Einstein-Rosen bridge (thus replacing the notion of "classical turning 
point" by that of "radius expectation value"). For the decreasing branch the 
same occurs on the "other" side of the Einstein- Rosen bridge. It is clear now 
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why the value of total mass-bare mass ratio in quantum case is greater than 
the corresponding classical value. The origin of this is just the replacing of 
the classical turning point by the mean value which is smaller thus giving 
rise to the increase in the total mass. 

Following such a line of reasoning we should conjecture that the maximal 
possible value of a total mass m for charge e and principal quantum number 
n fixed corresponds to the situation when the mean value of the radius of 
the shell lies at the event horizon making its collapse possible. The further 
increase in the bare mass m will lead to the wormhole with smaller total 
mass. 

Generalizing, we can introduce the notion of the quantum black hole 
states. For given values of electrical charge e and principal quantum number 
n the quantum black hole state is the state with maximal possible total mass 
m. Thus, for the quantum black hole states we have, instead of inequality 
( |155| ) the following equation 



3{nMy - AkM^c" + - An^ = 0, (157) 
the solution of which subject to the condition kM^ — > is 



kM^ = -e^ + -Ve^ + 12n2. (158) 
3 3 ^ ^ 



For the black hole mass spectrum we get eventually 
1 

rriBH 

'K 




(159) 



n = 1,2,3, 



For the uncharged black hole with e = we have immediately 

2 



m 



4 

27 



VnMpi. (160) 



Two limiting cases are of interest. In the first case 

<1, 



12722 

V 27 
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and the charge of the black hole gives rise only to the small corrections 
to the mass spectrum of the Schwarzschild black hole. The corresponding 
expressions for the bare mass M and the ratio m/M are, respectively, 



/4 V^/ V3e2 1 



m 




4 



M V 3 \ 12 n 96 n 

In the second case 



VSe^ 5 e' . /,„^x 

T7T--7^Z^ + ---]- (162) 



lei / n 



12n2 

2 



»l, 



The corresponding expressions for M and m/M are, respectively, 

M = 




M 2e4 ^ ^ 

It is interesting that in this second limiting case the black hole mass takes 
values nearly the mass of a classical extreme Reissner-Nordstrom black hole. 

It is easily seen that in both limiting cases | Am\/m <^ 1, where | A m| is 
the difference in masses for nearby energy (mass) levels. So, these limits are 
essentially quasiclassical ones, and the corresponding mass spectra should be 
compatible with the existence of Hawking's radiation. This is indeed the case 
as will be shown in the next section. 



42 



XII. Quantum Black Holes and Hawking's Radiation. 



In 1972 J.Bekenstein |^0[| recognized that the classical law that the area 
of a black hole horizon can not decrease has the striking resemblance of the 
second law of thermodynamics that the entropy of an isolated system can not 
decrease. There began the thermodynamical era of the black hole physics. 
The four laws of black hole mechanics has been formulated and proved [0 . 

This era became a quantum era in 1974 after S.W. Hawking discovered 
that due to quantum fluctuations black holes emit a blackbody radiation with 
the temperature which only by a numerical factor differs from Bekenstein's 
temperature, and the corresponding black hole entropy is just one fourth of 
the event horizon area. The starting point of our further investigation is 
the first law of the black hole thermodynamics which reads (for a spherically 
symmetric charged black hole) 

dm = ^QdA + (l>de, (165) 

where dm, dA, de are changes in a total mass m, event horizon dimension- 
less area A and electric charge e of the black hole, respectively, 6 is the 
temperature, $ is a Coulomb potential at the horizon, and 

4:71 r, , ^ 

A = —rl, 166 

K 

$ = (167) 
2 

e = — j-(r+ - Km), (168) 

r_|_ = Km + V K'^m'^ — Ke^. (169) 

The temperature of the Reissner-Nordstrom black hole is low in two different 
cases: when the mass of the black hole goes to infinity (for fixed charge) and 
when it goes to the minimal possible value y/Hm = e (extreme black hole). 
In both cases the back reaction of the radiation emitted is negligible, so the 
approximation of the static purely Reissner-Nordstrom black hole is valid. 
Moreover, it is just the quasiclassical limit of the (future) complete quantum 
picture in the sense that {\Am\/m<^ 1). The subsequent line of reasoning is 
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after V.Mukhanov []T2|. Here we extend his arguments to include the charge 
black holes, in particular, the nearly extreme Reissner- Nordstrom black holes. 

Let us suppose that the black holes are, in fact, quantizes, so their total 
mass and, hence, the area of the horizon are function of an integer number 
n. And let the black hole emits the quantum energy uj due to transition from 
the n-th state to the {n — l)-th state. The "typical" energy of this quantum 
minus the work done against the Coulomb attraction should be proportional 
to the temperature, i.e. 

ujn,n-i - ^de = TcaQ, (170) 

where de is the charge of the quantum in question, and a is a slowly varying 
function of the integer quantum number n. Inserting, then, Eqn. (|17CI|) into 
Eqn.( |165| ) we have {ujn,n-i = dm) 

dA = 4:na, (171) 



or 

dA = Airadn, (172) 

Integrating we get 

A = 47ran + AnC, (173) 

where C is a constant of integration, and a is some another slowly varying 
function of n. We consider two limits, n ^ C and n <^C . For the first case 
we assume a has the following expansion 

C C"^ 

a = ai+A- + 7i— = (174) 
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- < 1, 
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and for the latter one 
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a = a2 + /?27^ + 727^ = • • • , (175) 
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Let us write the expression for a black hole mass m in terms of a horizon 



area A making use of Eqns.(166) and {^6_ 



m 



1 t + 
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Inserting then Eqn.( PT^ ) for A, using expansions 
get 



(176) 

and (|176|) for a we 
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and 
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Comparing this with the expansion obtained from the black hole spectrum, 
we have 
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(179) 



Comparing values for ai and a2 we see that a is indeed a slowly varying 
function. 

Thus, we see that in the limiting case of low temperature our black hole 
spectrum is compatible with the existence of Hawking's radiation. Finally, 
we would like to note that our black hole spectrum obeys also the third law 
of thermodynamics. Indeed, for a nearly extreme Reissner-Nordstrom black 
hole the lowest energy level (for n = 1) always exceeds the critical value 
(= IcI/a/k), so the zero temperature state is not accessible. 
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Discussion. 



In the concluding remarks we would like to discuss shortly the obtained 
results. 

1. The appearance of the Schroedinger equation in finite differences is 
interesting by itself. And the fact that such an equation with Coulomb 
potential appeared to be exactly solvable seems miraculous. We believe that 
this is not accidental and there is some some symmetries that made the 
solution possible. Anyway, this point needs further investigation. 

2. The mass spectrum for the self-gravitating thin shells can be obtained 
from Dirac's spectrum for hydrogen atom by putting zero the so called radial 
quantum number (e.i. for the critical value of the total angular momentum). 

3. To extract the black hole spectrum from the thin shell spectrum we in- 
troduced the definition of the quantum black holes as a marginal case between 
quantum thin shells which do not collapse and wormholes. In quasiclassical 
regime this definition can be interpreted in that way that the expectation 
value of the radius of the shell is equal to the Schwartzschild radius. It is 
interesting that the special role played by the shells "lying" at the black hole 
horizon is also manifested in the processes of vacuum phase transitions which 
result in a black hole creation. In paper [|17| it was shown in a quasiclassi- 
cal approximation that the 0(3)-symmetrical vacuum phase transition with 
a black hole initially present or a blak hole formation from "nothing" may 
take place only if the new vacuum bubble shell is created exactly at the event 
horizon. 

4. We would like to emphasize once more the similarities an differences of 
our picture from the hydrogen atom models. In classical physics the electron 
in a hydrogen atom collapses radiating away electromagnetic energy. In the 
quantum mechanical picture it is stable unless something will force it to emit 
a photon and change its orbit, and there exists the lowest possible energy 
level which is by definition stable. Our self-gravitating shell collapses with- 
out any radiation in classical physics. The quantum mechanical picture is 
analogous to that of hydrogen atom. There exists also the lowest possible 
energy level. But unlike the hydrogen atom spectrum, our spectrum depends 
not only on the charge (which is already quantized for the electron) but also 
on a continuous parameter M, the bare mass of the shell. Moreover, we know 
that in the quasiclassical limit the black hole emits particles, thus acting as 
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a source of a blackbody radiation. This radiation is caused by the vacuum 
fluctuations of different fields and by existence of the event (apparent) hori- 
zon. It seems that the radiation is an inherent feature of the gravitational 
collapse. And to solve the problem of gravitational collapse we need to take 
into account the vacuum fluctuations and spontaneous emission caused by 
them. This will require the self-consistent solution of the gravitational and 
field theoretical (or string theory) problems. 

5. We found that each eigenvalue of a total mass (energy) is infinitely 
degenerate. Clearly, the origin of such a degeneracy lies in freezing all but 
one dynamical degrees of freedom. Our experience in quantum mechanics 
suggests that after restoring the frozen degrees of freedom this degeneracy 
disappears. As a result we expect that the "true" spectrum for thin shell 
masses will depend on more than one quantum number. Due to infinite 
degeneracy of the "frozen" spectrum the "true" one will, in fact, be quasi- 
continuous. And, as a consequence, the Hawking's evaporation spectrum will 
be quasicontinuous as well. The existence of the ground state with nonzero 
mass means that the evaporation process should stop when the mass of a 
black hole approaches the minimum possible value (~ 0.9Mpi) . And we 
may conjecture that the final state of the black hole evaporation will be 
characterized by the minimal value of mass but different wave functions of 
a ground state depending on details of mass distribution and evaporation 
process. We may call these wave functions quantum black hole hairs. This 
unexpected feature of the black hole mass spectrum (its infinite degeneracy) 
reminds the problem of so called hidden parameters in quantum mechanics 
and may be helpful in resolving the well known information paradox in a 
black hole physics. 

6. The very fact that our formulas may be obtained in the quasiclassical 
limits from the very general thermodynamical grounds encourages us to be- 
lieve that their range of validity is not confined to the specific model used 
for derivation. 7. There exists the minimal possible mass for the black hole. 
For the uncharged case it is slightly less that the Planckian mass. Such a 
minimal black hole mass was predicted by M.A.Markov |TB|. It is also the 



same for small electric charges and tends to the extreme value for charges 
comparable (in suitable units) to the black hole mass. For mass values less 
that this minimum there are no black hole states at all, and for all known 
elementary particles this is the case. It may appear that from the point of 
view of the spacetime structure all of them are on the "other" side of the 



47 



Einstein-Rosen bridge, forming wormholes. 
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